LOCAL EXTINCTION IN CONTINUOUS-STATE BRANCHING 
PROCESSES WITH IMMIGRATION 



CLEMENT FOUCART AND GERONIMO URIBE BRAVO 



Abstract. The purpose of this article is to observe that the zero sets of continuous state 
branching processes with immigration (CBI) arc infinitely divisible regenerative sets. Indeed, 
they can be constructed by the procedure of random cutouts introduced by Mandelbrot in 
1972. We then show how very precise information about the zero sets of CBI can be obtained 
in terms of the branching and immigrating mechanism. 

I. Introduction 

The problem of characterizing the zero set of a real-valued random process is, in general, 
not straightforward. When the process is a one- dimensional diffusion, several methods have 
been developed to study their zero sets. The problem is rather involved when dealing with 
Markov processes with jumps and this topic remains highly studied. We refer for instance 
to the recent survey of Xiao [Xia04] where fractal properties for Levy processes and other 
Markov processes are discussed. In this paper, we characterize the zero set of continuous-state 
branching processes with immigration (called also the CBI processes). 

Fundamental results on CBI processes, including their characterization as the large pop- 
ulation limit of Galton- Watson processes with immigration and the complete determination 
of the generator, are obtained by Kawazu and Watanabe [KW71]. Since then this class of 
processes have been studied extensively in several directions. 

The time evolution of the CBI process in general incorporates two kinds of dynamics; the one 
describing the reproduction and the other being due to immigration. Kawazu and Watanabe 
[KW71] show that the law of a CBI is characterized by the Laplace exponents ^ and $ 
of two independent Levy processes: a spectrally positive Levy process (which describes the 
reproduction) and a subordinator (which describes the immigration). 

One of consequences of introducing immigration is that zero is not an absorbing but a 
reflecting state. We provide necessary and sufficient conditions for zero to be polar, transient 
or recurrent. These results are obtained thanks to a connection between the zero set of a CBI 
and a random cutout set (see Mandelbrot [Man72]). Let Y be a CBI process started at zero 
associated to (\&, $). Define the random set 



Z := {t>0;Y t = 0}, 
and denote by v s the solution of the differential equation 

dv s T . . _ 

-j— = —W{v s ) and v 0+ = oo. 

Our main result is: 
Theorem 1. 

i) Z = {0} if and only if f^ exp [— J"" <f>(v s )ds\ du = oo. 
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ii) If Jq 1 exp [— J" &(v s )ds\ du < oo then 

a) T/ie random set Z is the closure of the range of a subordinator with Laplace ex- 
ponent 



L{q) 



e qt exp 



e 9 *exp 



§(v s )ds dt 



du I dt 



b) The random set Z has almost surely a positive Lebesgue measure if and only if 
Ie° wfsjds < oo (in that case, we say that the zero set is heavy, otherwise the set 
is light.) 

c) The random set Z is almost surely the union of closed non-empty intervals if and 
only if $ is the Laplace exponent of a compound Poisson process. 

Under an assumption of regular variation on the ratio R : u t— >■ we are able to 

give more details on the zero set. Loosely speaking, the index of regularity of the map R at 
+00 denoted by p measures the strengh of the immigration over the reproduction. We prove 
for instance that if p > — 1 then is polar. When p = — 1, some new constants are involved in 
the nature of the state zero. We get an upper and a lower bound for the Hausdorff dimension 
of the zero set. 

The connection made between random covering of the real line and the zero set of a CBI may 
be in some cases extended to another celebrated class of processes, the so-called generalized 
Ornstein-Uhlenbeck processes (OU processes). Besides obtaining results to generalized OU 
processes which are CBI processes, we provide a characterization of the zero set of Ornstein- 
Uhlenbeck processes driven by stable Levy processes. 



The paper is organized as follows. In Section 2, we recall the definition of the continuous 
branching processes, some fundamental properties and a key result on random covering. In 
Section 3, we study the zero set of a CBI started at 0. We state a corollary of Theorem 1 
and a general property about the fractal dimensions of the zero set. In Section 4, we establish 
precisely the connection with the random cutout set by applying the spine decomposition of 
a CBI (Theorem 9). In Section 5, we deal with a ratio function R with regular variation. 
Regular variation at +00 leads us to simple conditions on the upper and lower indices at +00 
of $ and ^ for polarity. Assumption of regular variation at 0+ leads us to conditions on the 
indices of $ and \1/ at for transience or recurrence. Finally in Section 6, we study the zero 
set of the Ornstein Uhlenbeck process driven by a stable Levy process. 



2. Preliminaries 

We begin by recalling some fundamental definitions and properties about continuous-state 
branching processes with immigration. We refer for instance the reader to Chapter 10 of 
Kyprianou [Kyp06] and Chapter 3 of Li [Lill] for a comprehensive survey on continuous-state 
branching processes. Then, we recall some theorems about random covering of the real line 
due to Mandelbrot [Man72], Fitzsimmons et al [FFS85]. The notion of random cutout set is 
at the core of the present work. 
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2.1. Continuous-state branching processes. Let x £ R+. A Markov process X(x) = 
(X t (x),t > 0), where X (x) = x is called a branching process in continuous time and contin- 
uous space (CB for short) if it satisfies the following property: For any y £ M + 

X(x + y)=X(x)+X(y) 

where X(y) is an independent copy of X(y). 

Let X{x) be a CB issued from x. There exists a unique triplet (d, a, v) with d, a > 0, and 
i/ a measure carried on M + satisfying 

DO 

(1 A x 2 )v(dx) < oo 

such that the Laplace exponent is given by 

E[e~ XXt{x) ] = exp(-OT t (A)) 

where the map t h-> ft (A) is the solution of the differential equation 

|f t (A) = -*(v t (A)), u (A) = A 

with 

2 poo 

(1) *(<?) = ^q 2 + dq+ (e- qx - 1 + gxl^i) 

The process is said to be critical, subcritical or supercritical according as ^'(0+) = 0, ^'(0+) > 
or \l/'(0+) < 0. For any x > 0, define the extinction time of the CB (X t (x),t > 0) by 

C := inf{t > 0;X t {x) = 0}. 

We have the following famous result (see for instance Theorem 3.5 and 3.8 pages 59-60 of 
[Lill]): 

P X [C < t] = exp (-xv t ) 
where v t = lim%(A). Recall the Grey's criterion: 

Theorem 2 (Grey, [Gre74]). The CB (X t ,t > 0) is absorbed in with positive probability if 
and only if there exists 9 > such that ^(z) > for z > 9 and 

dq 

< oo. 



Under that integrability condition, the real number v :=lim I v t £ [0, oof is the largest root of 

t—>oo 

the equation ^>(x) = and 

F x [( < oo] = exp (—xv). 

In the (sub) critical case v = and the CB(^>) is absorbed at zero almost surely. In the 
supercritical case, v > 0. 

Under the Grey's condition, the map (v 8 , s > 0) is the minimal solution of the differential 
equation (see Corollary 3.11 of [Lill]): 

(2) -^ = -9{v e ),Vo+ = oo. 
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2.2. Continuous-state branching processes with immigration. A continuous-state branch- 
ing processes with immigration (CBI for short) started at x is a Markov process Y(x) = 
(Y t (x),t > 0) satisfying the following property: for any y G K + 

Y(x + y) = Y(x)+X(y) 

where X(y) is an independent CB with mechanism Any CBI process is characterized by 
two functions of the variable q > 0: 

i r°° 

V(q) = dq+ -a 2 q 2 + J (e~ qu - 1 + qul uem )Mdu) 

POO 

$(q) = Pq+ / {l-e-« u )v (du) 
Jo 

where a 2 , (3 > and uq, v\ are two Levy measures such that 

/•oo roo 

l (1 A u)u (du) < oo and / (1 A u 2 )i>i(du) < oo. 
Jo Jo 

The measure V\ is the Levy measure of a spectrally positive Levy process which characterizes 

the reproduction. The measure z/ characterizes the jumps of the subordinator that describes 

the arrival of immigrants in the population. The non-negative constants a 2 and j3 correspond 

respectively to the continuous reproduction and the continuous immigration. Let F x be the 

law of a CBI (Y t (x),t > 0) started at x, and denote by E x the associated expectation. The 

law of (Yt(x),t > 0) can then be characterized by the Laplace transform of its marginal as 

follows: for every q > and x G M + , 



where 



^ x [e- qYt ] = exp (-xv t (q) - ^ $(v s (q))ds 

vt(q) =Q~ ®Mq))ds. 
Jo 



The pair $) is known as the branching-immigrating mechanism. A CBI process (Y tl t > 
0) is said to be conservative if for every t > and x G [0, oo[, f x [Y t < oo] = 1. A result of 
Kawazu and Watanabe [KW71] states that (Y t ,t > 0) is conservative if and only if for every 
e > 0, 

-dq = oo. 



o Mi)\ 

See Theorem 10.3 in Kyprianou [Kyp06] for a proof. 

Contrary to the simpler setting of continuous branching processes without immigration, 
non-degenerate stationary laws may appear. The following theorem provides a necessary and 
sufficient condition for the CBI to have a stationary distribution. Properties of the support of 
its stationary law have been studied by Keller- Ressel and Mijatovic in [KRM12]. 

Theorem 3 (Pinsky [Pin72], Li [Lill]). The CBI $) process has a stationary law if and 
only if 



tf'(0+) > and / -^-du < oo. 
Jo 



If the process is subcritical (*f>'(Q+) > 0), the convergence of this integral is equivalent to the 
following log- condition 

log(x)ui(dx) < oo. 



/ 

Jx>l 
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The main ingredient of the proofs provided in this work relies on a connection between the 
zero set of the CBI and a particular random set, called random cutout set. 



2.3. Random covering of the half real line. We recall here the definition of a random 
cutout set studied first by Mandelbrot [Man 72]. The notion of random cutout set is crucial 
in our study. For a lecture note presentation of the main theorems regarding random cutouts, 
we refer to the course of Bertoin [Ber99]. The main results we shall use in this paper may 
be found in. Consider a a-finite measure \x on R +) denote its tail fj,([x, oo[) by p,(x) for all 
x G E+. Let M be a Poisson point process on M + x E + with intensity dt <g> fi. Denote by 
(U,Xi)i e x hs atoms: 

We recover the half line by the intervals ]tj, tj + Xi[, i G Z. The set of uncovered point is 



1Z = [0,oo) -{J\ti,ti + Xi[. 



The measure \i is called the cutting measure. A random set 1Z is said to be a random cutout 
set if it is obtained by such a construction. Any random cutout set is a regenerative set (that 
is the closure of the range of a subordinator) . 

Theorem 4 (Fitzsimmons, Fristedt and Shepp, [FFS85]). 

If J exp (j p,(s)ds\ dt = oo then 1Z = {0} a.s. 

Otherwise 7Z is the closure of the image of a subordinator with Laplace exponent L given by 

= f e~ qt exp f / p,(s)ds J dt. 
HQ) Jo \Jt J 

Moreover if exp [ J 1 p,(s)dsj = oo then the set 1Z is unbounded. 



These random cutout sets have been intensively studied. As mentioned in the Introduction, 
we shall prove that the zero set of a CBI process is a random cutout set. Along the article, we 
shall use results on their geometry and refer the reader to the monography of Bertoin [Ber99] . 
For instance, we mention that Theorem 4 matches with Theorem 7.2 of [Ber99]. 

We end this section with the notion of infinitely divisibility for a regenerative set. As we 
shall see, this notion is underlying in all our work. A regenerative set TZ is said to be infinitely 
divisible if for any n > 1, there exist n independent identically distributed regenerative sets 
(7Z\ 1 < % < n) such that 



1Z 1 ={\1Z\ 

i=l 



Any random-cutout set is infinitely divisible. To the best of our knowledge, the converse is 
still an open question. 
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3. The zero set of continuous-state branching processes with immigration 

We shall focus on a CBI(\I/, <3>) started at 0. We assume that \l/ satisfies the Grey's condition. 
As we shall see in Section 4, this is of no relevance for the study of the zero set. By a slight 
abuse of notation, we denote by (Y t ,t > 0) the process (Y t (0),t > 0). Define the random set 



Z:= {t>0;Y t = 0}. 
Proposition 5. The random set Z is infinitely divisible. 

Proof. Let n > 1. Consider n independent CBI processes (Y l ; 1 < i < n) started at 
with branching mechanism \1/ and immigrating mechanism ^<3>. We have clearly the following 
equalities in law 



(Xut>0) = >0), 



and therefore 

n 

Z l = p|{t > 0;F/ = 0}. □ 
i=i 

If the set Z is not reduced to {0} and bounded (respectively unbounded), the state is 
transient (respectively recurrent). Recall that the state is said to be polar if 

F x [3t >0,Y t = 0] = 

for any i/0. 

Proposition 6. // the process is supercritical then is polar or transient. 

Proof. Assume that is not polar then by Theorem 1, we have 

-i 



exp 



®{v s )ds 



du = oo. 



To show that is transient, it suffices to show that the underlying subordinator is killed. This 
is equivalent to L(0) > 0, and therefore to 

exp I / &(v s )ds ) dt < oo. 



Since the process is supercritical (v t ,t > 0) converges as t — > oo toward a positive value v 
(the largest root of \l/). The map $ is continuous, then ($(t> t ),t > 1) also converges and is 
bounded. Denote by C a lower bound of &(v t ). We have f* &(v s )ds > C(t — 1) and then 

/oo / rl \ roo 

exp (j^ $>(v s )ds \ dt < J exp (-C(t - 1)) dt < oo. □ 

The next result is a corollary of Theorem 1. Since the proposition above deals with the 
supercritical case, we focus on (sub)-critical reproduction mechanism when studying the re- 
currence and transience of (statements ii) and iii) below). 

Corollary 7. Let 9 > such that ^(u) > for all u> 9. 

i) The state is polar if and only if 

1 



POO 




/ ex p 
le 


\I'TA 



Assume further that ^ is (sub) -critical, the state is 



dz = oo. 
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ii) transient if and only if 



poo 




/ ex p 
le 


\FTA 



9(z) 



dz < oo and 



iii) recurrent if and only if 



exp 



-dz < oo and 



exp 



exp 



g(f) 



r/.r 



< oo. 



oo. 



Remark 3.1. We may compare the integral conditions with those for the continuous branching 
process without immigration. The first statement has to be compared with the Grey's condition. 
The behavior of ^ in corresponds with the explosion of the CB (and also of the CBI). 
Assume that the process is conservative. This holds if and only if f = oo, if moreover 



t(x) 

Jo §(s)^ s < 00 ^ en ^ e secon d condition in statement iii) is always satisfied. Note that if the 
process is not conservative, obviously the state is either polar or transient. 

Proof. Proof of Corollary 7. Statements of the corollary are easily obtained by substitution in 
the integrals appearing in Theorem 1. Recall that vq+ = oo, the first statement is equivalent 
to the statement i) of Theorem 1 by considering u = v s . The transience or the recurrence 
of zero hold respectively if the set Z is bounded or unbounded. As already explained, this 
is equivalent for the underlying subordinator to be killed or not. If the state is not polar, 
then the zero set is bounded if and only if L(0) > 0. By the same substitution in the integral 

condition of Theorem 1, we get Q(v t )dt = JJ 1 exp 

the limit of (v t ,t > 0) when t goes to oo and equals to since we focus on the (sub)critical 
case. The constant v\ does not play any role here and may be replaced by 9. □ 



Jx w(s) 



dx 



The constant v is 



Recall the definition of the lower and upper box-counting dimension. See for instance 
Section 3 of [Xia04] or Chapter 5 of [Ber99]. For every non-empty bounded subset E of 1R + , 
let N e (E) be the smallest number of intervals of length e needed to cover E. The upper and 
lower box-counting dimension of E are defined as 

hg(N e (E)) 



and 



dim(_E) := limsup- 

log(l/e) 



dim(E):=liminf^™. 

e^o log(l/e) 



We provide now our last general result on the zero set. The proof is postponed at the end of 
Section (4). 

Lemma 8. The random set Z has the following upper and lower box-counting dimension. 
For every t > 0, 

1 



and 



dim(Zn [0,t]) 



dim(Z n [0,t]) 



lim inf 

u^O log( 



limsup - — 

u^O l0g( 



-ds a.s. 



-ds a.s. 



Moreover if Z is bounded almost surely then the law of the last zero of (Y t ,t > 0), 

goo := sup{s > 0; s G Z} 
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is given by 




i 




Pf^oo e dt] = k 1 exp 



&(v s )ds ] dt — k 1 exp 



du dt, 



with k the renormalization constant. 

Remark 3.2. The lower box-counting dimension and the Hausdorff dimension (denoted by 
dimn) of a regenerative set coincides almost surely (see Corollary 5.3 of [Ber99]/ 

4. Analysis of the zero set: spine decomposition and random covering 

The objective of this section is to recall a construction of CBI processes which will then be 
used to prove Theorem 1. 

A CB(\&) process reaches zero with positive probability if and only if the branching mech- 
anism \& satisfies the Grey's condition. Intuitively, a CBI(\I/,$) cannot touch zero, unless 
the corresponding CB(\I/) can reach zero. We begin by establishing rigorously this idea. Let 
x > 0; the branching property stated in Subsection (2.2) provides that the CBI^, $) process 
Y(x) satisfies 



If x > and Y(x) reaches zero with positive probability then plainly the CB(\1/) X(x) reaches 
also zero with positive probability. Theorem 2 ensures then that \1/ verifies Grey's condition. 
In order to handle the case x = 0, we note that $ ^ implies that v t (X) > and so Y t (0) > 
with positive probability for any t > 0. If Y(0) reaches zero with positive probability then 
it does so after time t for some t > 0. Applying the Markov property at time t on the set 
Y t (0) > 0, we are reduced to the previous case and conclude that Grey's condition holds. Thus, 
we see that imposing Grey's condition is of no relevance to study zero sets of CBI processes. 

In order to prove Theorem 1, we now recall a Poissonian decomposition of CBI processes, 
also called the spine decomposition, which is obvious in the simpler setting of Galton- Watson 
processes in immigration and indeed has been considered by a number of authors in the 
continuous setting. For example, [PY82] use it to obtain a decomposition of Bessel bridges or 
in [Eva93] it allows representations of sup erpro cesses conditioned on non-extinction, a work 
which has continued in [EW03, FL04]. 

The spine decomposition of CBI processes is found in subsection 2.1 of [CR11] or as a 
particular case of the construction of subsection 2.2 of [AD09]. It is based on the N-measures 
constructed for superprocesses in [DK04] and specialized to the case of CBI processes in 
Theorem 1.1 of [CR11]. We now give a streamlined exposition of the construction of our 
specialized N- measure, based on [PY82], by making the assumption 



Assumption (3) is exactly Grey's criterion for reachability of zero for CB(\1/) processes. 

Starting in [PY82], a Poisson process representation for (continuous) CBI processes has 
been achieved by means of a a-finite measure which can be understood as the excursion law of 
a CB(\1/) although, in general, one is not able to concatenate excursions to obtain a recurrent 
extension of a CB(\E'). See for instance Subsection 3.1 of [Lam02]. 

For a CB(^) we have 



Y(x) = Y(0)+X(x). 



(3) 




= e 



xv t (\) 
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for any x > 0. Hence, A i— >■ v t (X) is the Laplace exponent of a subordinator; Under (3), it is a 
driftless subordinator since 

lim^ = 0. 

A— >oo A 

See for instance Corollary 3.11 page 61 of [Lill]. (This is the fundamental simplification in 
comparison with [CR11].) Let rj t stand for the Levy measure of the Laplace exponent v t so 
that 



v t 



(A) = / (1 - e~ Xx ) r) t (dx) 
Jo 



by the composition property v t o v s (X) = v t+s (X) we can see that r\ t is an entrance law for the 
semigroup of the CB(\I/) so that we may consider the a-finite measure Q on the space of cadlag 
excursions starting and ending at zero characterized by its finite-dimensional marginals: for 
< ti < ■ ■ • < t n 

Q(X h E dX!,...,X tn e dx n ) = l^.^oTjt^dX!) P^_ tl (Xi,dx 2 ) ■ ■ ■ P^t^iXn^^Xn) . 

This can be thought of as an excursion law for the CB(\I / ). This is because as t — > 0, r\ t — > 5q 
since v t (X) — > X so that under Q, X = 0. Also, under Q, X is Markovian and with the 
semigroup of the CB(^). To characterize the image of the length of the excursions under Q, 
say (, notice that 

Q(C >t)= lim Q(( > t + h) 

h—i-0+ 

= lim f F*((>t) Vh (dx) 

h->0+ J 

= lim f (1 - e~ xvt ) rj h (dx) 

h-^0+ J v ) i \ i 

= lim v h (v t ) = v t . 

h->0+ 

Now use $ and Q to construct the cx-finite measure N on excursion space by means of 

N{dX) = 0Q{dX) + J v {dx)Fl{dX) 
which will be the intensity of the Poisson point process 

Finally, let 

Theorem 9. Y is a CBI(^, $) which starts at 0. 

This is what we refer to as a spine decomposition of a CBI process. 
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Proof. We start with a preliminary computation which uses the exponential formula for Poisson 
point processes: 

E(e- xy t)=E (exp-A £ X t O 

V ieX;ti<t J 

= exp (- j N(l - e" AXt - s )^^ 

= exp jT* J™ (i _ e -^-*W) ^(rfy) d s -p jf y (1 _ e -^) ^(dar) 
= exp ^-^ $(u t _ s (A)) ds^ = exp ^- ^ $(u 8 (A)) ds^ . 



Hence, at least Y has the correct one- dimensional distributions. To prove that Z is a CBI, 
we need to compute conditional expectations. Let 



Tt — o 



Then 



E 



exp | ^ w 4+ 

t<*i<t+S 



E 

exp 
exp 



\ t<ti<t+S J 



t+s 



<&{v t+s _ r (\)) dr 



v u (\) du . 



On the other hand: 



E 



exp 



t+S—ti 



exp(-5>i^ s (A)) 
exp(-v s (A)F 4 ) . 



Hence 



E(e 



-XY t - 



\T t ) = exp ( -Y t v s (X) 



<5>(v u (\))du 



□ 



Remark 4.1. The Grey's condition fulfilled here by the map ^ is not necessary for this spine 
decomposition to hold. We refer the reader to Theorem 1.1 of [CR11]. 

Proof of Theorem 1. From this spine decomposition we establish a useful connection between 
the random covering procedure and the set of the zeros of the CBI. Results due to Fitzsimmons 
et al. [FFS85], recalled in Subsection (2.3), will allow us to study the set Z. Denote for each 
i G X, Q = inf{t > 0; X\ = 0}. By standard properties of random Poisson measures, the 
random measure Yli^i ^(u,d) is a Poisson random measure with intensity dt <g> N(£ G dt). 
Plainly, we have 



{t>0,Y t = 0} 



\ 



U [ti,ti + ti[ U |J + 



ieX;X^>0 
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The following crucial equality holds 

(4) Z = R + \U ie i]t h t i + { i [. 

Notice that some U may belong to M + \ U ie i]ti,ti + Q[ and not to {t > 0, Y t = 0}. However 
such point ti is a left accumulation point of the set {t > 0, Y t = 0}. Namely, consider ti £ Z, 
assume by contradiction that there exists e > such that [ti — e, tj[n{i > 0; Y t = 0} = 0, then 
taking the closure we have [t, — e, ti] HZ = ® which is impossible since £ Hence, equality 
(4) holds almost surely. Moreover, we have 

POO 

N(C > t) = /3Q[C >t]+ ¥ x (( > t)u (dx) 



(3v t + 



1 - exp(-x^)) u (dx) 



In order to establish the statement, we can directly use the results of Subsection (2.3). Part 
i) and statement a) of ii) readily follow from Theorem 4, taking 

m = n(c >t) = $00 

for all t > 0. We get Z = {0} if and only if 





/"U 








du = oo 


/o 







Using Equation ((2)), the substitution z = v t and the condition w + = oo yields 

1 











/ exp 


-J §(v s )ds 


du = exp 


Jv 1 . 


Jo 




J VI 



dz = oo. 



Statement b) follows from Proposition 1 of [FFS85] (see also Corollary 7.3 of [Ber99]). Namely 
the set Z has a positive Lebesgue measure if and only if 

$(v s )ds < oo. 



(s A l)n{ds) 
By the same substitution, this yields 



ds < oo. 



The last statement d) follows from Corollary 2 of [FFS85]. Namely, the necessary and sufficient 
condition for the uncovered set to be a union of intervals is that the intensity measure N(£ £ dt) 
has a finite mass. Therefore, assume that $(u t ) - — > c < oo. We have v t — > oo and 



${vt)=Pv t + (1 



-v t x 



By monotone convergence 



)u (dx). 



M dx ) ^o([0,oo[). 



Thus we get the conditions (3 = and uq([0, oo[) < oo. 



□ 



We now provide the proof of Lemma 8. The arguments are directly those used for the fractal 
dimensions of random cutout set. 
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Proof of Lemma 8. Combining Theorem 5.1 and Corollary 7.6 of [Ber99], we get that for every 
t > 0, 



dim(Z n [0,t]) = 1 - liminf h ®&*) ds as 



log(^) 
and 

r ,> [ ®(v s )ds 

dim(Z n 0,t ) = 1 - limsup Ju \ ' .{ a.s. 

u^O log(-) 

We get the statement by substitution u = v s . The second statement concerning the largest 
zero is a direct application of Corollary 7.4 of [Ber99]. □ 

5. Regularly varying branching-immigrating mechanisms 

In this section, we shall focus on specific mechanisms \& and $ such that the function ratio 
R(x) = ^||y is regularly varying at oo. The index of R, called in the sequel, p may be 
interpreted as representing the strength of the immigration over the reproduction. For sake of 
concisness, we only work with a critical branching mechanism A remarkable phenomenon 
occurs when p = — 1. In such case, other quantities are involved. Surprisingly the quantities 
r := liminf sR(s) and f := limsup sR(s) play a crucial role. 

5.1. Polarity of zero and Hausdorff dimension of the zero set. We first recall the 
definition of the upper and lower indices at infinity of a Laplace exponent (we refer the reader 
to the seminal work of Blumenthal and Getoor [BG61] and to Duquesne and Le Gall's work 
[DLG05] page 557). 

log(^(X)) 

Ind(fr) = supia > 0; lim ^(\)\- a = oo} = liminf f V ,\" , 

A-s-oo A->oo log(A) 

IHd(tf) = inf{a > 0; lim *(A)A~ Q = 0} = limsup lo f ^ffl . 

A->oo A ^oo log(A) 

This definition holds also for the Laplace exponent of a subordinator $, replace only \1/ by $ 
above. 

Lemma 10. Assume that the function x i— >■ R(x) = |^ is regularly varying at oo with index 
p, then p = Ind (Q) — Ind(^f). Furthermore, we have Ind (<&) — Ind(^f) = Ind(<&) — Ind j^f). 

Proof. We will prove that 



kgf^rr) — ► (lnd($)-IndW)log(A). 

\ R{X) J x->oo 



First write 

,'R(Xx) \ log($(x)) - log(*(x)) 

log 



R(x) J log(x) 



x J \x J \og(R(x)) 



Recall that since R is regularly varying with index p, there is a slowly varying function I such 
that R(x) = x p l(x). Therefore 



I IV log (-) lQg(jR(Ax)) - log 
k x J \x J \og(R(x)) \x 



log(/(a:)//(Ax))-plog(A) 
p\og(x) + \og(l(x)) 
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Since log ( 777^ ) — > 0, the limit of the quantity in the right hand side of the last equality is 
the same as 

plog(A) log(x) log(A) 

plog(x) +log(/(x)) ~~ 1+ ll£g0(£»" 

P log(x) 

Since the function I is slowly varying, we have that — ► (see for instance Proposition 

1.3.6-(i) of [BGT89]). Finally, since 



R(x) 

we get p = Ind($) - Tnd(\^) = Tnd($) - Ind(W). □ 

Theorem 11. i) If Ind (&) > Ind{^>) — 1, then is polar. 

ii) If Ind (&) < Ind(ty) — 1, then is not polar and the zero set is heavy. 

iii) If Ind (<&) = Ind{^) — 1. Define the quantities 

r := hmsup s - , . 

and 

r := hmmt s . 

Soi/i quantities belong to (0, 00]. 

a) J/r < Ind (ty) — 1, then is not polar. Moreover ifr > then the zero set is light. 

b) If r> Indfif) — 1 then is polar. 

Proof. By assumption the map R is regularly varying at +00 with index p, assume that 
p > — 1, therefore z H- ff R{u)du is regularly varying with index p + 1 and by Karamata's 
theorem (see for instance Proposition 1.5.8 of [BGT89]) we have 

(5) For all z >9, [ R{u)du = C + z p+1 l(z) 

Je 

with C a constant and I a slowly varying function. By Corollary 7, we have to study the 
following integral 

f°° ( r \ 1 

I := exp I / R{u)du dz. 



i) If p > —1, by Proposition 1.3.6-iii) of [BGT89], we have z p+l l(z) — > +00, moreover 

z— >oo 

observing that ^f(z) = 0(z 2 ), we have clearly that X = 00. 

ii) If p < —1, by Proposition 1.5.10 of [BGT89], we have J e °° R{u)du < 00, then X < 00 
and furthermore by statement ii)-b of Theorem 1, the zero set has a positive Lebesgue 
measure. 

iii) In the case p = — 1, the map R is regularly varying with index —1. We mention that 
integrals of such regularly varying functions yields to de Haan functions (see Chapter 
3 of [BGT89]). 

(a) Let e > 0. For s large enough, R(s) < - + - . Therefore for large enough z, 

f r \ 1 z 7+e 

exp / R(s)ds — - < C— T . 
\Je J * (2) V(z) 
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Moreover V(z) > CV^W -6 and thus 

Therefore, as e is arbitrary small for large values of z, 

z r-]nd(*)+2e dz < x if F _ Ind(^r) < _ L 

This implies that is not polar. Assume r > for s large enough, R(s) > r/2s, 
we have clearly J e °° R(s)ds = oo. 

(b) For s large enough, sR(s) > r — e and then exp ( J g 2 R(s)ds) > Cfj^y By the 
same reasoning, we have ^(z) < Q' z Ind ^ +e and we get 



Since e is arbitrary small, J e °° z r 2e d2 = oo if r — Ind( l I/) > — 1. □ 

Remark 5.1. It should be highlighted that we may study the polarity of zero without the regular 
variation assumption. For instance, we easily get sufficient conditions for zero to be polar or 
not. Namely if Ind (§) > Ind(^>) — 1 then is polar, and if Ind(§) < Ind (^) — 1 then is 
not polar. However the case Ind(&) > Ind (ty) — 1 and Ind (&) < Ind(^>) — 1 may occur and is 
rather involved. Our assumption of regular variation avoids such cases. 

Proposition 12. Assume that f < Ind (^) — 1 then for all t > 0, 



1 - = < dim(Z n [0, t]) < 1 - =- 



Ind(^) - 1 ~ Ind{^) 
and _ 

< dimjZ n [0,t]) < 1 



Indm - 1 v Indm - 1 ' 

Before tackling to the proof, we need a lemma. 
Lemma 13 (Duquesne Le Gall, Lemma 5.6 of [DLG05]). The following equalities hold 

liminf j = =1 

log(l/t) Ind(^>) - 1 

and 

log(ut) 1 

hmsup- — - . . = - — — — -. 

t^o log(l/t) IndW-l 

Proof. Duquesne and Le Gall established 



liminf , < 



't->o log(l/t) ~ IndO) - 1 
and 

log(vt) 1 

llHlSUp- . . < - — — — -. 

t _,o log(l/t) " Ind(^) - 1 
The other inequalities are also true. We provide here a proof. Denote r] := Ind(\l/), there exist 
e > arbitrarly small such that 

q^ +e) ^(q) — ► 0. 

q—^oo 
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For all C > 0, for large enough q we have > -^7. Therefore for small enough t, we have 



dq + f°° dq C 1-m-e 



This yields 



We then have 



( - J < log (j-^ — - ) 



liminf- — l ) s > 



't->o log(l/t) 77 — 1 + e" 

As e is arbitrarly small, the result follows and the first equality in the statement is established. 
We study now the second equality. Denote 7 := Ind(\l/). Let 7' > 7. There exists a sequence 
(u n ,n > 1) such that u n — > 00 and ty(2u n ) < 2 7 u} L . For u e]0,2tt n ], we have by convexity 



*(«) 1 1 (2u„) an( ^ therefore ty(u) < 2 7 ' 1 mm 7 ' 1 when u g]0,2m„]. Define the function F( 



a 



f°° -^rrr, we have 



Then 



u mi nf lo f/fW)<i im lo f/f'"»»<y-i. 

a-s-oo log(a) n->oo log(w n ) 



As page 592 of [DLG05], observe that by definition of v t : 

logK) \ ^ . f log(l/F(q)) 

hmsup- — tttt =hmmt — . 

t^o log(l/t) / 1^00 fo#(a) 

We deduce that 

log(u t ) 1 
hmsup- — -— - > - — -, 
t^ log(l/t) 7 - 1 

by letting 7' go to 7, we obtain the wished inequality. □ 
Proof of Proposition 12. Let e > 0, by assumption, we have for C large enough 

sup sR(s) < r + e and inf sR(s) > r — e. 

se[c,oo[ se[c,oo[ 

Therefore, 

(r-e)[log(^)-log(C)]< jf ^t4^< (r + 6)[logK)-log(C)]. 
By using the previous lemma and Lemma 8, we plainly get 

1 r + e — < dmT(z n [0, t}) < 1 



Ind(tf) - 1 ~ Ind(tf) - 1 

As e is arbitrary small, we get the statement. Same arguments hold for the lower box-counting 
dimension. □ 
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5.2. Recurrence and regular variation at 0+. To study the recurrence of zero, we need 
information on the behaviour of the map R in the neighbourhood of 0+ (see statement iii) in 
Corollary 7). In the same vein as our previous study of the polarity, a natural assumption is 
to consider the map R with regular variation at 0+. 

In order to state our result, we need to introduce the lower and upper indices of a Laplace 
exponent at 0+: 

ind(tf) = supia > 0; lim^(A)A" a = 0} = liminf f l ,\' J 

A->0 A->0 log(A) 

logWA)) 

ind(tf) = inf{a > 0; lim^(A)A^ Q = oo} = limsup , ,\* . 

a^o A ^ log(A) 

Remark 5.2. If J 0+ -j^y = oo then we can easily observe that 1 < mrf (f) < ind{9) < 2. 

Moreover if the function \I/ is regularly varying at 0+ then ind(^>) = ind (^). Note that a non 
trivial example for which ind (^>) = 1 is given by ty(x) = xlog(x) (however this example does 
not satisfy the Grey's condition). 

Exactly with the same arguments as in Lemma 10, we have the following property. 

Lemma 14. Assume that R is regularly varying at 0+ with index k, then we have: 

k = ind(<f>) - ind(^) = ind(§) - ind{^). 

Theorem 15. If the map R is regularly varying at +oo with index p and at 0+ with index k, 
then if p < — 1; 

i) is transient if k < —1 

ii) is recurrent if k > — 1 . 

If p = — 1 and r < Ind (fy) — 1 then is recurrent if k > —1, transient if k < —1. 
Proof. We assume that is not polar, then we have to study the following integral 

dx 



J ■= exp 



R(u)du 



y(x) 



Assume that R is regularly varying at with index k G R, then by an easy adaptation of 
Propositions 1.5.8 and 1.5.10 of [BGT89] to the setting of regular variation at 0+, we get 

lim f 9 R(u)du — \ °° H K < \ 

x^0 Jx K ' \ < OO if K > — 1. 

More precisely, if k < —1, we have J 9 R(u)du = x K+l l{x) with / a slowly varying function at 
0+. 

• Assume n < —1. Let e > 0, we have for small enough x, ^(x) > Cx md ^ +e . We shall 
prove that the map 

exp(— x K+1 l(x)) 
x i — y = 

j,ind(*)+e 

is bounded in the neigbourhood of 0. This will provide us that J is finite. Taking the 
logarithm, we have 

/ exp(-x K+1 l(x))\ K+ll/ , / ,— — , v . 1 log(x) 

log —, : = -x K+1 l(x) 1 + (md ^) + e)— - ,/ 

aJnd *)+« J \ K x ) x 
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On the one hand, x K+1 l(x) — > +00, on the other hand the map x 1— > lo ,^ x ? is slowly 



varying at 0+, and thus using Potter's bound, we have x — > 0. Finally 



x^O 



and the map is bounded. We deduce that J < 00 and then is transient 
if k > — 1, provided that J 0+ = 00, we have J = 00 and is recurrent. 
We deal now with the 1. Denote 

7t := limsup xR(x) > 0. 

x^O 

Let e > 0, for u small enough, we have 

R{u) < jj + I and then - f x R(u)du > -(k + e) log(x) + C 
Observing that 

exp - / R{u)du ) > C^—j- — > 



*(^) \ Jx ) ' ^( X ) " ^( X ) 

we get 

Finally, if k = — 1 (and is not polar) then is recurrent. □ 

We study in the sequel the specific case of stable and gamma mechanisms. 

Remark 5.3. If k > — 1 then J 0+ R(u)du < 00 and by Theorem 3, the CBI has a stationary 
law. 

5.2.1. Stable and gamma mechanisms. Consider of the form ^>(q) = dq a and = 

d'q 13 with a G (1,2], /3 G (0,1] and d,d' G (0, 00). The CBI process (Y t ,t > 0) associated 
is said to be stable. Obviously, we have Ind(\I/) = a, Ind($) = 0. The map R is regularly 
varying at 0+ and at +00 with index p = n = (3 — a. We can thus apply the previous results 
in Theorem 11 and Theorem 15. 

• If (3 > a — 1, then is polar, 

• if f3 < a — 1, then is transient and the zero set is heavy, 

• if /? = a — 1, we have r := lim xR(x) = 4" and k — — 1. Two cases may occur. 

x— >oo 

— If > ct — 1 then is polar, 

— if 4 < a — 1 then is recurrent and by Proposition 12 

dim H (Zn[0,t]) = l-^—^-. 

a — 1 d 

Notice that in this stable framework, we cannot have p < — 1 and k > — 1. 

Proposition 16. Let (Y t ,t > 0) denotes a stable critical CBI started at with parameters 
a, (3 satisfying (3 = a — 1. 



Z = {t>0;Y t = 0} = {a tl t>0} 



with (a t ,t > 0) a stable subordinator with index 7 = 1 — 



1 d' 

q-1 d 
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Proof. In order to get that the subordinator (at,t > 0) is a 7-stable one, we shall use self- 
similarity result. Indeed Patie showed in [Pat09] that this is the only self-similar CBI process. 
The self-similarity followed by inspection. Namely, we have 



E[e- gYt ] = exp (- J <S>(v s (q))ds 



with v s (q) = q[l + d(l — a)q a H] a - 1 . An easy computation yields that the processes (kY t , t > 
0) and (Yjfea-it, t > 0) have the same law. We deduce that the regenerative set Z is self-similar 

meaning that for any k > 0, kZ l = Z. The only regenerative sets satisfying this property 
are the closure of the range of stable subordinator (see for instance Section 3.1.1 of [Ber99]). 
Proposition 12 provides the Hausdorff dimension of the set and therefore the index of stability 
(see Theorem 5.1 of [Ber99]). □ 

Consider the immigration mechanism $(0) = r T/£^, 9 l ~ ^r^- The subordinator with 

Laplace exponent $ is called Lamperti stable subordinator (see for instance [CPP10]). The 
map $ is regularly varying at +00 with index /3 and at 0+ with index 1. Assume that the 
reproduction mechanism \I/ is a-stable, as previously we can apply Theorem 11 and 15. We 
easily get 

• if P > a — 1 then is polar, 

• if /3 < a — 1 then is recurrent and the zero set is heavy, 

• if j3 — a — 1, we have r := lim sR(s) = dr ^_^ and k = 1 — a > — 1. Two cases may 
occur. 

— if d < ft"t then is polar 

— r(a) 1 

— if d > prHr then is recurrent. 

r(a) 

Finally, by Proposition 12 we get 

1 



dim H (Z n [0,t]) 



T(a)d 

Last, consider now the case of a Gamma immigration mechanism. Let a > 0, b > be the 
parameters and 

POO 

$(x) = a log(l + x/b) = (1 - e~ xu )au- x e~ hu du. 







We may observe that <3> is slowly varying at 00 and regularly varying at 0+ with index 1. In 
this framework, considering a stable reproduction: ty(x) = x a , we get: R : x 1— > |^jy regularly 
varying at with index k = 1 — a and at +00 with index p = —a. Applying Theorems 11 and 
15 we have that if a G (1,2], then k > — 1 and is recurrent. Moreover the zero set is heavy. 



6. Ornstein-Uhlenbeck processes 

We recall here basics on Markov processes of Ornstein-Uhlenbeck type. Contrary to CBI 
processes, these processes in general take values in R. Consider x G M, 7 G 1R + , and (A t ,t > 0) 
a one-dimensional Levy process with characteristic function given by 77 such that 

E[e izAt ] = exp{t V {z)), 



a 



2 



r){z) = - — z l + ibz+ / (e lzx - 1 - izxl { \ x \< l} )v{dx) 
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where a > 0, b G K and / K (|a;| 2 A l)r/(da;) < oo. A process (X t , t > 0) is said to be an OU type 
process if it satisfies the following equation 

X t =x-7 f X s ds + A t . 
Jo 

Generalized Ornstein-Uhlenbeck processes valued in R + belong actually to the class of the CBI 
processes. When A is subordinator with Levy measure v and drift d, the Ornstein-Uhlenbeck is 
a CBI with = 72; and $(-2) = cte + /g 00 ^ — e~ zx )v(dx). In this case, the Grey's condition 
is not fulfilled and clearly the state is polar. Most of the Ornstein-Uhlenbeck processes 
are not CBI processes (we refer for instance to the discussion about their stationary laws in 
Proposition 4.7 in Keller-Ressel and Mijatovic [KRM12]). However, there is an interesting 
class of OU processes whose zero sets are random cutout sets: a-stable Ornstein-Uhlenbeck 
processes. 

Theorem 17. Let (X t ,t > 0) be an a-stable Ornstein-Uhlenbeck process started at 0. If 
a G (0, 1] then the zero set of X equals {0} almost surely. If a & (1, 2] then the zero set of 
X is a random cutout set whose cutting measure has density with respect to Lebesgue measure 
given by 

z^ (1 - (3)e z /(e z - ifdz 

where = 1 — 1/a. 

Proof. Consider the process 

Xt = e ^A e \ta. 

Because A is self-similar, it follows that X is stationary. Using integration by parts, it follows 
that 

dX t = e~ xt dA e x ta - \tX t dt. 
Then, by the fact that A is a stable Levy process, we infer that J Q e~ Xs dA e \ as has the same 
law as A. Hence X is the stationary version of the Ornstein-Uhlenbeck process driven by A. 
On the other hand, the zero set of X is then the logarithm of the zero set of A, which is a 
self-similar regenerative set and therefore a random cutout set. Actually, the zero set of A is 
empty if a G (0, 1] (cf. [Ber96, page 63]) while if a G (1, 2] then it has the law of the (closure) 
of the image of a /3-stable subordinator with j3 — 1 — 1/a (see [BG62], [Get63] and the latter's 
note added in proof). 
Let 

iez 

be a Poisson point process with intensity dt® where fi(dx) = (1 — f3)/x 2 dx. Then the zero 
set of A has the law of a random cutout based on H and so the zero set of X is the random 
cutout (on M.) obtained by removing the intervals 

(Si, Si + Zi) = (log(ti), l0g(*i + Xi)). 

Namely, we have 

Z:={teR;X t = 0} = R- \J]s i: Si + z\ 
The intensity of the point process 

is the image of the measure dt ® n(dx) by (t,x) 1— > (s,z) = (log(t), log(£ + x)). A notable 
cancellation occurs and it is found to be equal to (1 — f3)e z /(e z — l) 2 dsdz. Using the identity 
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established between the cutting measures of random cutouts on (—00, 00) and random cutouts 
on (0,oo) (cf. Theorem 2 of [FFS85]) we see that the zero set of the a-stable Ornstein- 
Uhlenbeck process (started at zero) is a random cutout set whose cutting measure has density 

(1 - f3)e z /{e z - ifdz. □ 

The cutting measure with density (1 — f3)e z /(e z — l) 2 was studied in detail in example 8 of 
[FFS85]. There, the authors show that the associated subordinators have zero drift and Levy 
measure v given by 



z/((x,oo)) 



c 



Note that the density of v can also be written as 

Ce x C'e^ x 
u( - dx ^ = ( e * _ dx = ( sinh (x/2)) 2 -/ 3 dX - 

Hence, in the special case of the Ornstein-Uhlenbeck process associated to Brownian motion, 
we recover the results of [PY97] (which go back to [HT91]). However, we also deduce that 
in the general stable case, the zero set is the image of the Lamperti stable subordinators 
introduced in our special case in [DMY07] and studied in general in [CPP10]. 



Corollary 18. If a G (1,2], the random set Z = {t > 0,X t = 0} is infinitely divisible. More- 
over Z is almost surely not bounded (0 is recurrent) and we have for all t > 

dim H (Zn [0,t]) = 1/a. 

Proof. We only have to give a proof for the Hausdorff dimension. The Laplace exponent of 
the Lamperti stable subordinator involved in Theorem 17 is k(A) = r(i-p)r(\) ( see Equation 
(27) in [FFS85]). Therefore dim H (Z n [0, t)) = Ind(re) = 1 - (3 = 1/a. □ 

Acknowledgment. CF would like to express his gratitude to Jean Bertoin and thanks 
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